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We have measured the isotope shift of the narrow quadrupole- allowed 5 2 5i/2 - 4 2 D 5 / 2 transition 
in 86 Sr + relative to the most abundant isotope 88 Sr + . This was accomplished using high-resolution 
laser spectroscopy of individual trapped ions, and the measured shift is Af^f c 'f s 6 = 570.281(4) MHz. 
We have also tested a recently developed and successful method for ab-initio calculation of isotope 
shifts in alkali-like atomic systems against this measurement, and our initial result of Ai/ 88 ^ 86 = 
457 (28) MHz is also presented. To our knowledge, this is the first high precision measurement 
and calculation of that isotope shift. While the measurement and the calculation are in broad 
agreement, there is a clear discrepancy between them, and we believe that the specific mass shift 
was underestimated in our calculation. Our measurement provides a stringent test for further 
refinements of theoretical isotope shift calculation methods for atomic systems with a single valence 
electron. 



I. INTRODUCTION 



Isotope shift measurements provide precise tests of our 
understanding of atomic and nuclear structure [1-6]. Ad- 
ditionally, precise knowledge of isotope shifts is useful in 
studying alternative explanations for observations of an 
apparent temporal variation of the fine structure con- 
stant [7]. The problem of isotope shift calculation is 
rich in relativistic many-body physics, and it has his- 
torically presented a strong challenge when considering 
atomic species with many electrons. Motivated by the 
aforementioned applications, there has been much recent 
work toward improving theoretical methods for calculat- 
ing isotope shifts [7-9]. Along with others [1, 5, 10], the 
measurement reported here provides a stringent test case 
for our calculation method [7] and other suitable methods 
in the context of ions with alkali-like electronic structure. 
In this paper, we present a high precision measurement 
of the A = 88 to A = 86 isotope shift in the S\/2 - -D5/2 
quadrupole transition of Sr II, a medium- weight ion, and 
we also test the ab-initio calculation method described 
in [7] with this measurement. The isotope shift of this 
quadrupole transition in Sr II has proven to be in a par- 
ticularly interesting regime for theoretical work, and, to 
our knowledge, the calculation presented here is a first 
reported calculation of this shift. 
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The experimental technique used here involves 
quantum-shelving spectroscopy of individual ions, in 
which kilohertz-level determinations are possible due to 
the trapping of an atom in a potential with a strength 
significantly larger than the recoil energy of the photons 
used to interrogate it [5]. The very high spectroscopic 
precision of the measurement is also afforded, in part, by 
the use of a dipole-forbidden transition with a linewidth 
that is approximately eight orders of magnitude smaller 
than the more commonly studied dipole-allowed transi- 
tions in Sr II [2, 3, 11-13]. With the natural linewidth 
of this transition being a mere 0.46 Hz [14, 15], the pre- 
cision of the present measurement is limited by the fre- 
quency stability of the probe laser used to interogate the 
quadrupole transition. To our knowledge, the most pre- 
cise experimental determination of the 88 Sr + - 86 Sr + iso- 
tope shift in the S\/ 2 - -D5/2 transition prior to the present 
one was made at the megahertz level [10] (though from 
the format in which the result is presented, this prior de- 
termination may only be at the 10 MHz level; no uncer- 
tainty is quoted). Our measurement improves upon this 
precision by more than two orders of magnitude. We also 
note that a similar measurement of 247.99(4) MHz for the 
88 Sr + - 87 Sr + isotope shift of the same quadrupole tran- 
sition has also been reported in the context of a 87 Sr + 
ion optical frequency standard [1]. 

Aside from providing a very stringent test for atomic 
and nuclear structure calculations in atomic systems pos- 
sessing alkali-like electronic structure, our measurement 
also provides a useful tag for identification of 86 Sr + in ion 
frequency standards or quantum information processing 
applications, especially when combined with absolute fre- 
quency measurements [16, 17] for the same transition in 
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FIG. 1: (Color online) Level diagram for even isotopes of 
singly ionized strontium (i.e., Sr II) with no net nuclear spin, 
88 Sr + and 86 Sr + in this context. The 6*1/2 level is the ground 
state of the single valence electron in these ions. The P levels 
decay rapidly, having lifetimes of approximately ten nanosec- 
onds, and the D levels are metastable (lifetime approximately 
half a second). The double-arrows indicate laser radiation 
used in the determination of the isotope shift in the 445 THz 
(red) quadrupole transition. The thick lines represent dipole- 
allowed (El) transitions while the thin line indicates the much 
weaker dipole-forbidden but quadrupole- allowed (E2) transi- 
tion of interest. The natural lifetimes r of the excited states 
are taken from [15]. The eigenstates of the optical qubit in 
the quantum information processing apparatus used to per- 
form this measurement are |0) = IS1/2) and |1) = [-D5/2}- 



the more abundant and commonly used isotope 88 Sr + . 
One application in which this tag might be useful is in 
investigations of sympathetic sideband cooling of one or 
more motional modes of a multi-isotope coulomb crystal, 
where one ion species is utilized for "memory qubits" and 
the other is used as "refrigerator" ions [18-23]. This may 
be of intereset as the isotope shift reported here may be 
large enough for use of laser radiation derived from the 
same 445 THz source to simultaneously perform both 
sympathetic near-ground-state cooling and coherent op- 
erations on memory qubits in a manner similar to the 
work of [22], i.e., with an optical qubit encoded in the 
S1/2 and -D5/2 levels (cf. Fig. 1). 

A detailed account of our experimental methods is pre- 
sented in section II followed by the results in section III. 
The discussion proceeds with pertinent details of the cal- 
culation of the quadrupole transition isotope shift in sec- 
tion IV. This report concludes with a comparison of and 
remarks on our theoretical and experimental results in 
section V. 



II. EXPERIMENTAL METHODS 

We trap Sr + ions in a four-rod radio frequency (RF) 
Paul trap, described previously [23, 24], with axial and 



radial trapping frequencies of ~ 500 kHz and ~ 1.2 MHz, 
respectively. Ions are loaded one-by-one from a neutral 
strontium atomic vapor via photoionization (PI) using a 
two step process to an autoionizing level with laser ra- 
diation at 650 THz and 740 THz [25]. Ions are Doppler 
cooled to approximately the Doppler limit using light of 
frequency 711 THz, slightly red-detuned from the 6\ / 2 - 
A/2 dipole-allowed cycling transition. We have noticed 
a small amount of isotope selectivity during ion loading 
by tuning the first-step PI laser slightly to the red by a 
few hundred megahertz so that we can load 86 Sr + from 
our unenriched strontium sample somewhat more often 
than would be expected from the natural abundances of 
strontium isotopes (9.86% 86 Sr+; 82.6% 88 Sr+). We oc- 
casionally load 87 Sr + as well. However, our current ap- 
paratus does not allow direct cooling or detection of this 
species due to its ground-state hyperfine splitting. The 
P1/2 level decays to the D 3 / 2 level, which has a 0.4 s 
lifetime, approximately 1 in 14 times. Consequently, a 
"repumping" laser at 274 THz continuously illuminates 
the ion to depopulate this level and return the valence 
electron to the cooling cycle (cf. Fig. 1). The ion is in 
the Lamb-Dicke regime [26] for the quadrupole-transition 
probe light at 445 THz, thus allowing high resolution 
first-order Doppler-shift-free spectroscopy on this nar- 
row transition via quantum shelving spectroscopy of the 
ion [27]. 

The 445 THz light is generated using a homebuilt 
external-cavity diode laser that is narrowed and stabi- 
lized with a Zerodur low-thermal-expansion optical ref- 
erence cavity. The cavity is housed in a vacuum chamber 
in a thermally controlled enclosure with temperature sta- 
bility at the millikelvin level, and the laser is locked to 
this cavity using the Pound-Drever-Hall technique [28]. 
Feedback is provided both to the external laser cavity 
grating and to the laser diode current to narrow and 
frequency lock the laser. This has allowed a demon- 
strated lower bound on the optical qubit (cf. Fig. 1) 
coherence time of 2.5(0.2) ms when using simple refo- 
cusing techniques to circumvent technical and environ- 
mental noise (e.g., magnetic field fluctuations at the ion 
location) [23] . On the much longer timescale of the sev- 
eral seconds that are required to perform precision spec- 
troscopy on this optical qubit, the probe laser frequency 
varies due to drifts in both the length of its reference 
cavity and the locking electronics. On this timescale, 
the demonstrated linewidth is still less than 1 kHz, thus 
providing for the observation of narrow quadrupole tran- 
sition spectral lines [23] . After each interrogation on this 
transition and subsequent state determination, the ion 
state is reset to the ground state Si/ 2 level by driving 
any population in the -D5/2 level to the P3/2 level from 
where it predominantly decays to the ground state. This 
"quenching" is accomplished using a 290 THz laser (cf. 
Fig. 1). By applying this radiation simultaneously with 
the 274 THz repumping light, the ion will be optically 
pumped to the S'1/2 state with very high probability in 
several tens of nanoseconds. The 445 THz spectroscopy 
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beam is pulsed using an acousto-optic modulator (AOM) 
in a double-pass configuration, which allows for scanning 
of its frequency over a few tens of megahertz. The cooling 
and quenching laser beams are switched off using shutters 
in order to avoid coupling of the D5/2 level to other lev- 
els, which would otherwise broaden and shift the narrow 
quadrupole transition while it is being probed. 

Each quantum shelving spectroscopy trial proceeds as 
follows. The ion is first Doppler cooled for approximately 
10 ms with the Doppler, repumping, and quenching laser 
radiation on. The Doppler radiation and the quencher 
are then shuttered off. Approximately 1 ms later, the 
narrow transition is probed via a pulse from the spec- 
troscopy laser of varying frequency detuning and typi- 
cally of 1 to 10 ms in duration. Another few to several 
milliseconds later, the Doppler laser is switched back on, 
and after less than 1 ms, fluorescence is detected using 
a photomultiplier tube for 5 ms with near unity state 
detection efficiency. If the ion was not excited by the 
spectroscopy pulse, it will fluoresce and a mean of ap- 
proximately 35 Poisson-distributed photon counts are ac- 
quired in this detection time. If the spectroscopy pulse 
succeeded in exciting the ion to the -D5/2 level where it 
then remains shelved (i.e., removed from the Dopplcr- 
cooling cycle) during detection, a mean of approximately 
1.8 photons is collected. This background is due pri- 
marily to 711 THz Doppler cooling light scattering from 
the trap electrodes into the fluorescence collection op- 
tics. Finally, the quencher beam is switched on to pump 
the ion back to the ground state in case it was shelved, 
and the ion is then re-cooled for a few milliseconds prior 
to commencement of the following trial. The repumper 
beam remains on during the entire experiment; the AC 
Stark shift and broadening due to this radiation on the 
narrow transition of interest are negligible because of its 
~ 16 THz detuning from the closest transition involving 
the -D5/2 level. 

At each frequency detuning of the probe laser, 100 tri- 
als were carried out. A photon count threshold, typically 
of 9 counts, was set in order to minimize error in distin- 
guishing dark signals from bright ones due to the overlap 
of their respective Poisson distributions [29], where tri- 
als with a photon count at or below the threshold were 
considered as dark and all others were bright. The ratio 
of trials where the ion was dark to the total number of 
trials at a given detuning gave the excitation probabil- 
ity at each frequency step. Spectra like that shown in 
Fig. 2 were constructed in this manner, and they display 
the Zeeman components of the S'1/2 - -D5/2 transition. 
The laser detuning that corresponded to the center fre- 
quency of a given quadrupole transition spectrum was 
determined as a parameter of a fit to the spectral data. 
The total magnetic field at the location of the ion was also 
determined from the Zeeman spectra, and each trial was 
triggered on the AC line to minimize effects from mag- 
netic field fluctuations at the line frequency of 60 Hz [36] . 
We carried out these determinations for both isotopes of 
interest, and comparison of the RF frequencies delivered 
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FIG. 2: (Color online) Sample spectrum of the S'1/2 - -D5/2 
transition in Sr II acquired using the quantum shelving tech- 
nique. The blue curve is a fit to eight Zeeman components of 
this quadrupole transition with the four most central Zeeman 
lines shown here. The dashed red line indicates the measured 
RF frequency delivered to the AOM at the transition center, 
which was determined as a parameter of the fit. 



to the AOM at the transition centers enabled the deter- 
mination of the isotope shift. 

The total shift is on the order of the approximately 
600 MHz free spectral range (FSR) of the 445 THz probe 
laser reference cavity, and the two different cavity fringes 
used to lock that laser for probing the two isotopes were 
in adjacent FSR's as a result. The full determination 
of the isotope shift therefore required a measurement of 
the FSR of the probe laser reference cavity as well. The 
total shift in frequency detuning is the AOM component 
of the shift Av^qm measured from the shelving method 
just described with subsequent subtraction of one FSR, 
Aupsr, of the reference cavity (i.e., the lock point for 
86 Sr + is approximately one cavity FSR red of that for 
88 Sr+): 
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(1) 



We note that the convention Av A ' A — v A — v A is fol- 
lowed here, where v A and v A are the energies of the 
transition of interest in the two isotopes with mass num- 
bers A' and A with A' > A. 

To determine Avfsr with some degree of precision, we 
used an electro-optic modulator (EOM) in an auxiliary 
measurement to impose frequency sidebands at approxi- 
mately half of the cavity FSR on the portion of the probe 
laser light that is coupled into the fixed-length reference 
cavity. While scanning the grating of the laser, we mon- 
itored the transmission through the cavity; by adjust- 
ing the EOM drive frequency such that the first-order 
blue sideband from one FSR overlapped the first-order 
red sideband from an adjacent FSR, we could determine 
Avfsr to a precision limited by the laser linewidth and 
the cavity finesse at 445 THz. We measured the reference 
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FIG. 3: (Color online) Measures of the AOM component 
Ai/'q* of the overall isotope shift. Black circles are mea- 
surements of the S — D transition in 86 Sr + , and red squares 
are measurements of the same transition in 88 Sr + . The AOM 
frequency is the RF frequency shift delivered to the probe 
light at the resonance line center not including the offset due 
to the free spectral range of the reference cavity that was used 
to lock and narrow the probe laser. The laser was locked one 
cavity FSR to the red when measuring the 86 Sr + transition 
compared to when measurements of the 88 Sr + transition were 
performed. The insets show the same points with zoomed 
ranges to clearly show the effect of cavity drift (note the dif- 
ferent ordinate scales). An overall frequency offset of 200 
MHz has been subtracted from all RF frequency values for 
readability. 



cavity FSR to be 

A^ F sfl = 599.301(2) MHz. 



(2) 



The isotope shift, as measured using the shelving 
method, is affected by an isothermal drift in the refer- 
ence cavity length. Through measurement of the tran- 
sition resonance in individual ions over several months, 
we have determined the long-term average of the cavity 
drift to be 52(1) kHz/d. On shorter time scales, how- 
ever, the cavity exhibits nonlinear drift, and it has been 
observed to drift at rates up to approximately six times 
the long-term drift rate over a few hours. To date, our 



isotope-shift measurements have been performed on only 
one ion at a time. This leads to the requirement that 
^ v avm mus t be determined from the difference of two 
measurements taken at different times, with the isotope 
shift determined as a parameter of a fit to the nonlinear 
cavity drift as explained below. Lack of detailed knowl- 
edge about this nonlinear drift as a function of time dom- 
inates the imprecision of our current determination. It 
should be noted that this drift corresponds to a maximum 
of approximately 3.6 Hz/s. As the frequency scan dura- 
tions were 2 to 4 minutes, we estimate the error from cav- 
ity drift during a scan to be a few hundred hertz at most, 
which is significantly below our uncertainty in center fre- 
quency determination at the few kilohertz level, and it 
has been ignored in this analysis. The cavity drift also 
affects Avpsr through the change in the length of the 
cavity, but the fractional length change is approximately 
1.2 x 10 -10 per day, leading to a nominal drift in the 
FSR of 3 mHz/h (up to ~ 9 mHz/h), which is well below 
our uncertainties during all measurements contributing 
to the present isotope shift determination. 

The data used in the analysis consist of (i) three mea- 
surements of the Si/ 2 - -D5/2 transition frequency in 
86 Sr + , (ii) three measurements of the same transition fre- 
quency in 88 Sr + , and (iii) the measurement of the probe 
laser reference cavity FSR, as described above. The six 
S1/2 - D5/2 measurements were obtained over a span 
of approximately four hours, a timescale over which the 
nonlinear reference cavity drift has a minimal effect on 
our uncertainty. The RF frequencies delivered to the 
probe light via the AOM at the centers of the acquired 
Zeeman spectra are plotted in Fig. 3 versus the relative 
time of each measurement (with an arbitrary zero of time 
at 00:00 local time on the day that the data were col- 
lected). The zoomed insets show the effect of cavity drift 
between measurements. 

To best account for the cavity drift when determining 
the AOM shift, Az/^, 8 ^ was introduced as a fit parame- 
ter, and a simple drift function with a discontinuity was 
fit to the six measurement points versus time. As the ex- 
act character of the cavity drift during data acquisition 
was not known, two different drift functions, one linear 
and one quadratic, were chosen: 



h{t) = (mit + 6i)e(t -*) + (mi* + 6i + Ai/^)e(t-to) (3) 
h(t) = (<Z2< 2 + m 2 t + b 2 )e(t Q - t) + (q 2 t 2 +TO 2 i + fo 2 + A^ 88 o 8 ^)O(i-7; ). (4) 

I 

Here O is the Heaviside unit step function, and the qi, not required [37]. 

rrii, and bi are fit parameters for function fi of time t. 

The discontinuity time to can be any time between the 

86 Sr + and 88 Sr + measurements, and it is chosen here to 

be 17 h. The reference cavity drift is sufficiently slow 

such that a higher than quadratic order fit function is 
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FIG. 4: (Color online) Fits to determine AOM component of 
isotope shift. Black circles are measurements of the S — D 
transition in 86 Sr + , and red squares are measurements of the 
same transition in 88 Sr + . The blue dashed line is a lin- 
ear fit yielding an AOM shift of 29.020(2) MHz, and the 
green solid line is a quadratic fit, yielding an AOM shift of 
29.019(2) MHz. The two panels are offset in the ordinate co- 
ordinate by 29.020 MHz, but are scaled identically (refer to 
left and right scales), as are the abscissas. An overall fre- 
quency offset of 200 MHz has been subtracted from all RF 
frequency values for readability. 



III. DISCUSSION OF EXPERIMENTAL 
RESULTS 

The linear fit to f l yields Avf^ = 29.0204(21) MHz 
with a reduced \ 2 (goodness-of-fit parameter) of 1.54 
(3 degrees of freedom). The quadratic fit to f 2 yields 



Ai/ 



88.86 
AOM 



29.0191(18) MHz with a reduced \ of 0.88 (2 
degrees of freedom); see Fig. 4. As the actual drift func- 
tion during these measurements is unknown, we weight 
the average of these two measurements of Av^qm °y ^ ne 
inverses of their respective goodness-of-fit parameters. 
This gives a value of Az/f^ = 29.020(3) MHz, where 
the uncertainty has been increased from 0.002 MHz to 
0.003 MHz to conservatively account for lack of knowl- 
edge of the drift function. With the value for Avp sr m 
equation (2), we calculate a value (cf. equation 1) for the 
quadrupole transition isotope shift to be 



A^ 8 n 8 oa 8 s 6 = 570.281(4)MHz. 



(5) 



To our knowledge, this is the first measurement of this 
shift to high precision, at a level of better than 1 part in 
10 5 , with the uncertainty dominated by nonlinear drift in 
the probe laser reference cavity. From equation (5) and 
the known value of the absolute frequency of the 5 2 Si/2 - 
4 2 D 5 /2 transition in 88 Sr + [17], we also infer an absolute 
frequency for this transition in 86 Sr + of 



v m = 444, 778, 473, 814(4)kHz. 



(6) 



IV. CALCULATION OF THE ISOTOPE SHIFT 

To an excellent approximation, isotope shifts in atomic 
transition frequencies arise from two sources [30-32]: (i) 
the mass shift (MS) and (ii) the field shift (FS). The MS 
is the difference in recoil energy between isotopes due to 
their differing nuclear masses; i.e., when conserving mo- 
mentum in the atomic systems, the eigenenergies of their 
electronic states shift to account for the differing recoil 
of the nuclei. The MS is itself comprised of the sum 
of two different effects referred to as the "normal mass 
shift" (NMS) and the "specific mass shift" (SMS). The 
NMS is due to a reduced mass correction which accounts 
for the finite mass of a real atomic nucleus compared to 
the infinitely massive nucleus of an idealized atom, and 
it is readily calculated. The SMS accounts for the mod- 
ification of the nuclear recoil energy due to correlations 
in the motions of all possible pairs of electrons orbiting 
the nucleus via conservation of momentum applied to the 
atomic system. Accurate calculation of this shift is a dif- 
ficult problem that depends heavily on the accuracy with 
which the multi-electron wavefunction for the atom can 
be determined. The FS is a consequence of the differ- 
ence in the neutron-number-dependent size and shape of 
the nuclear electric charge distribution between isotopes; 
i.e., the eigenenergies of the electronic states adjust to the 
unique electric potentials of their nuclei. Measurement 
of the FS requires the use of an electronic transition that 
involves a significant change in the electronic wavefunc- 
tion density that penetrates the nuclear volume (or in 
approximate terms, A \ip(r = 0)| =/= 0, where r = is 
the nuclear center and ip is the electronic wavefunction). 
The electronic states most sensitive to nuclear screening 
are the S states followed by the P states, with states in 
all other angular momentum levels having negligible pen- 
etration of the nuclear volume [30]. Thus, in principle, 
the quadrupole transition utilized here should be a sen- 
sitive probe for measuring the FS in this medium-weight 
atomic ion where neither the SMS or FS may be treated 
as negligible a priori. 

The difference in the transition frequency Av A ' A be- 
tween isotopes with mass numbers A' and A can be ex- 
pressed as [7] 



&v A '- A = (k NM s + ksMs) 



-F8(r 2 ) A ' A . (7) 



Here, the normal mass shift constant is 



WMS — 



1823' 



where v is the transition frequency of the idealized atom, 
1823 refers to the ratio of the atomic mass unit to the 



A' A 



electron mass, and <5(r 2 \" is the relative difference in 
the mean squared nuclear charge radius between the two 
nuclear species. The SMS and FS constants, ksMS an d 
F, must be calculated. 

The Sr II ion considered in this paper can be treated 
as a single valence electron above a closed-shell electronic 
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core. Therefore the methods of ref. [7] for computing F 
and ksMS are appropriate. Full details are given in [7]; 
here we only present an outline. The general idea is to 
include an additional term which represents the SMS or 
FS — a "finite field" — in an energy calculation from the 
very beginning. We rescale the finite field using a param- 
eter A, chosen to make the effect large but numerically 
stable. The energy calculation is then repeated for differ- 
ent values of A and the numerical constant is extracted 
as the gradient, i.e. ksMS = dE/dX. 

Calculation of the energy is done with the AMBiT 
code [8, 9]. We use the Dirac-Fock method to generate 
the core orbitals and their potential. Valence and vir- 
tual orbitals are then constructed via a large i?-spline 
basis [33] so that core-valence correlations can be in- 
cluded using many-body perturbation theory (MBPT). 
We calculate these correlations to second order in the 
operator H — Hjjf-, the difference between the exact and 
Dirac-Fock Hamiltonians. Higher-order correlations are 
estimated using Brucckncr orbitals. An operator £ is cre- 
ated such that the second-order MBPT correction to a 
valence orbital \n) is n|£|n). Then E is added 

to the exchange potential in the Dirac-Fock Hamiltonian 
for the valence electron and this new potential is used 
to create the single-electron Brueckner orbital. We use 
the difference between the second-order solution and the 
Brueckner orbital solution as an estimate of the errors 
associated with higher-order correlations. 

For the Si m - -D5/2 transition in Sr 11 we obtain 

fc/VMS = —244 GHz • amu, 
ksMS = -If 49 (50) GHz • amu, and 
F = -4616 (272) MHz/fm 2 . 

From this intermediate result, we calculate a total mass 
shift of 368 (13) MHz and a field shift of 89 (f 5) MHz, 
where we use J(r 2 ) 88,86 = -0.0549 (17) fm 2 from [34]. 
Adding these two contributions (cf. equation 7), we fi- 
nally obtain 

A^f= 457(28)MHz. (8) 

Since we have measurements for Az/^f e 'f,f (this work) 
and A^ 88 c 'fJ [f] for the same quadrupole transition in 
Sr II, we may also use equation (7) and the appropriate 

SOr*) ' x from [34] to determine ksms ari d F experimen- 
tally. These constants are determined to be 

ksMS = —1660 (20) GHz • amu and 
F = -1200(100) MHz/fm 2 . 

We see that our ab-initio calculation has significantly un- 
derestimated the magnitude of the specific mass shift in 
this ion, a discrepancy which is not accounted for by the 
error estimation method. 



V. SUMMARY AND CONCLUDING REMARKS 

We have presented a high precision measurement, to 
better than 1 part in 10 5 , of the 88 Sr+ - 86 Sr+ isotope 
shift of the 5 2 S , i/ 2 - 4 2 D 5/ / 2 quadrupole transition of Si- 
ll that is consistent with the result reported in [10] and 
which improves upon that reported precision by more 
than two orders of magnitude. The recently developed 
and successful method [7] for the ab-initio calculation of 
isotope shifts in single electron atomic systems has also 
been tested against our new measurement, and we use the 
present measurement in conjunction with an analogous 
measurement of the 88 Sr + - 87 Sr + isotope shift [f ] to per- 
form an experimental check on our calculation. Sr II is a 
medium-weight atomic ion that presents a somewhat dif- 
ficult case for theoretical calculations of its isotope shifts 
since the SMS and FS are both large, and neither one 
may be regarded as negligible, thus making it difficult 
to separate these two quantities in equation (7) without 
some independent information about either one of them. 
There is a clear discrepancy between our calculated iso- 
tope shift and the measured values. The error is most 
likely to be in the SMS, since the operator for this effect, 
^i<jPi • Pj; Pfc being the momentum of electron k, is very 
sensitive to details of the electronic wavefunctions. Rel- 
ativistic corrections to the isotope shift are at the level 
of a few percent in light atoms [35], and should not be 
much larger in our case. Therefore, the most likely reason 
for the discrepancy is higher-order core-valence correla- 
tions (beyond second order) that are not included in our 
calculation. Interestingly, our method for estimating the 
magnitude of higher-order correlations (using Brueckner 
orbitals) has failed in this case: the calculated ksms is 
31% too small, but our predicted error was 6%. Never- 
theless, the calculation is in broad agreement with the 
measurement, which is a promising initial result when 
one considers that the specific mass shift is unusually 
dominant in this case; it is approximately 7 times the 
size of the normal mass shift. Our high precision mea- 
surement presents a strong test for further refinements in 
theoretical iosotope shift calculation techniques for alkali- 
like atomic systems. 
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